Thermally- Assisted Spin- Transfer Torque Magnetization Reversal in Uniaxial 

Nanomagnets 
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We simulate the stochastic Landau-Lifshitz-Gilbert (LLG) dynamics of a uniaxial nanomagnet 
out to sub-millisecond timescales using a graphical processing unit based micromagnetic code and 
determine the effect of geometrical tilts between the spin-current and uniaxial anisotropy axes on 
the thermally assisted reversal dynamics. The asymptotic behavior of the switching time (J — ► 0, 
(r) oc exp(— £(1 — /) 2 )) is approached gradually, indicating a broad crossover regime between ballistic 
and thermally assisted spin transfer reversal. Interestingly, the mean switching time is shown to 
be nearly independent of the angle between the spin current and magnet's uniaxial axes. These 
results have important implications for modeling the energetics of thermally assisted magnetization 
reversal of spin transfer magnetic random access memory bit cells. 



A spin-polarized current passing through a small mag- 
netic conductor is known to deposit its spin-angular mo- 
mentum into the magnetic system [l|, Q ■ This in turn 
causes the magnetic moment to precess and in some cases 
even switch direction. This has led to sweeping advances 
in the field of spintronics through the development and 
study of spin-valves and magnetic tunnel junctions (see, 
for example, HI). A thorough understanding of these 
phenomena, however, requires taking into account the 
effect of thermal fluctuations inducing diffusion in the 
moment of a magnetic material. This is of particular ex- 
perimental relevance since spin-transfer effects on nano- 
magnets are often conducted at low currents, where noise 
is expected to dominate. Recent debate in the literature 
over the proper exponential scaling behavior between 
mean switching time and current elucidates how the ther- 
mally assisted properties of even the simplest magnetic 
setups leave much to be understood 043 • By simulating 
reversal dynamics of a uniaxial macrospin model sub- 
ject to spin-torque, we demonstrate how harnessing the 
vast computational parallelization capabilities intrinsic in 
Graphics Processing Unit (GPU) technology for numeri- 
cal modeling can allow a deeper probing of the thermally 
activated regime. This leads us to a further understand- 
ing of the energetics at play in small magnetic systems 
subject to spin-torque effects. 

The standard theoretical approach used to study the 
dynamics of small magnetic structures has been to treat 
the magnetic system as a macrospin in the spirit of 
Brown [9( , thus modeling the dynamics as those of a mag- 
netic moment undergoing damped precession under the 
action of an effective field. After Slonczewski one 
assumes that the magnetic body absorbs spin-angular 
momentum perpendicularly to its orientation. As such, 
spin-torque effects are generally taken into account phe- 
nomenologically by modifying the macrospin's Landau- 
Lifshitz-Gilbert (LLG) dynamical equation with the ad- 



dition of the spin-torque term. The LLG equation for a 
then reads: 



rh = — 7'm x H c ft — aj'm x (m x H e ff) 

— 7'sm x (m x rip) + 7'Qsm x n p , (1) 

where 7' = 7/(1 + a 2 ), 7 is the gyromagnetic ratio, a 
the Landau damping coefficient and s = {fi/2e)rjj is the 
spin-angular momentum deposited per unit time with 
V = (^t — Jl)/(Jf + Jl)> the spin-polarization factor of 
incident current J. The last two terms describe a vector 
torque generated by current polarized in the direction n p . 
The effective field H c ff is immediately obtained once the 
energy landscape of the magnetic sample is specified. In 
the simplified case of a uniaxial monodomain magnet, it 
is characterized by the projection of the magnetization 
onto the uniaxial anistropy axis: 



U(m) = -K(m ■ h K f m H. 



(2) 
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with K = (1/2)MsHkV the uniaxial anisotropy of a 
magnet of volume V and anisotropy field Hk, saturation 
magnetization Mg, applied field H, and hx the unit vec- 
tor representing the orientation of the uniaxial anisotropy 
(easy) axis. 

Of special interest is the interplay between spin-torque 
and thermal noise effects. A spin current also induces 
magnetization fluctuations. But for typical experimen- 
tal situations this noise is far smaller than the thermal 
noise, which we discuss below. Starting from (1), ther- 
mal fluctuations are modeled by considering the effect 
of a random applied magnetic field with mean zero and 
variance enforced by the fluctuation dissipation theorem. 
Doing this leads to a stochastic LLG equation with mul- 
tiplicative noise 0, [T^] • Experimentally, one is generally 
interested in the magnetization reversal properties of a 
thin film magnet. For currents less than a threshold cur- 
rent (to be defined below), the switching in such magnetic 
systems appears to be dominated by thermal-activation 
related reversal processes. Fitting experimental data to 
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theory is complicated by several factors. First, the va- 
lidity of the macrospin model is questionable for samples 
typically studied (see, for example, HU), which are so 
far larger than the exchange length of transition metals 
magnets (~ 5 nm). Second, even within the macrospin 
approximation it is not clear what the relevant energetics 
are for thermally assisted spin transfer torque reversal. 
This is a consequence of the inherently non-conservative 
nature of the added spin-torque term which does not al- 
low the construction of a generalized energy landscape 
suitable for Kramer's theory analysis. 

The energetics of thermally assisted reversal in the 
presence of spin transfer torques for a macrospin is the fo- 
cus our analysis. Here, theoretical progress has been hin- 
dered by the computational power needed to run numer- 
ical simulations to the desired degree of accuracy. The 
LLG equation, modified as a set of coupled stochastic 
equations can be studied in one of two ways: either by 
concentrating on its associated Fokker-Planck (FP) equa- 
tion or by constructing a stochastic Langevin integrator 
to gather sufficient statistics on reversal. In either case 
simulations have thus far been unable to extrapolate to 
long enough times to approach the thermally activated 
regime. If on one hand, it is hard to maintain conver- 
gence of direct numerical solutions of the FP equation 
out to long timescales, gathering sufficient statistics to 
make proper use of a Langevin integrator has proven to 
be just as difficult a task. A recent paper by Taniguchi 
and Imamura Q, suggests that the proper scaling be- 
havior between mean switching time and current should 
be log((r)) oc —(1 — J) 2 , with I = s/{oiHk) the normal- 
ized current. This contradicts previous work proposing a 
proportionality to — (1 — I) Nonetheless, numerical 

simulations have not yet been able to settle this disagree- 
ment. 

To elucidate the most basic properties of thermally- 
assisted magnetization reversal under spin-transfer 
torque, we choose to consider a magnetic system solely 
characterized by its uniaxial anisotropy and focus on its 
thermally assisted dynamics in the absence of any applied 
fields. The energy landscape is given by the first term on 
the right hand side of (2). The equilibrium states in this 
landscape are simply given by parallel and anti-parallel 
configurations of the magnetization along the uniaxial 
anisotropy axis. We consider our magnetic sample to be 
initially magnetized in thermal equilibrium with one of its 
two possible stable states (i.e. m = +n#- or m = —hie)- 
Upon turning on a (positive) current, spin-torque effects 
will attempt to drive the magnetization toward align- 
ing with the polarization axis rip. If the current is then 
switched off again, the magnetization will relax toward 
the easy axis once more. To capture this behavior, it is 
very helpful to write down the dynamical equation for 
the projection q = m • % of the magnetization along the 
easy axis. For temperatures such that £ = K/kT 3> I, 
upon turning off the current, the sign of q will specify to 
a high degree of accuracy in which state the magnetiza- 
tion will relax. We call this the "projectional dynamics" , 
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FIG. 1: Mean switching time versus current in the sub-critical 
low current regime (/ < 1). Times are shown in units of (s-T) 
where T stands for Tesla: real time is obtained upon division 
by Hk- The red and green line are obtained by fitting the 
long time data to the functional form (r) = C exp(— £(1 — J)' 1 ), 
where C is fitted numerically and fi is an exponent, either I 
or 2. The red curve (/i = 2) can be seen to approximate the 
numerical data asymptotically better than the green curve 
(fx = I). The theoretical curve (in light blue) represents 
the asymptotic switching scaling given by Eqn. 4. The dis- 
crepancy between numerical data and theoretical asymptotic 
scaling of the mean switching time seems to indicate a broad 
crossover regime between ballistic and thermally assisted dy- 



which read: 



a 



{n z I + q)(l- q 2 ) + n z Iq(q - — ) 

n z 



a 



+ aIn x m y + x j—{l-q 2 )oW. 



(3) 



For definitiness, we have conventionally chosen the z-axis 
to be collinear with the polarizer (n v = z) and the easy 
axis to be lying in the z-x plane. To simplify the nota- 
tion, then, n x and n z are simply the components of fif. 
Furthermore, a natural time unit t = j'Hxt has been 
introduced. The last term appearing in the expression 
models the effects of thermal noise whose stochastic con- 
tribution is intended according to Stratonovich calculus 
(indicated by the symbol o) and W is a standard mean 
0, variance 1 Wiener process 12 1. 

In the simplified collinear case, where polarizer and 
uniaxial axes are aligned, the projectional dynamics are 
the same as the regular dynamics for the m z component 
of the magnetization, and are decoupled from both m x 



and 



The problem is effectively reduced to that of 



a ID stochastic differential equation characterizing the 
switching behavior over all possible currents. It is easily 
seen that 1=1 behaves as a crossover value above which 
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all dynamics prefer the parallel state independently of 
the noise, and this value of I is hereafter denoted as the 
"critical current". Analogously, the thermal regime is 
characterized by sub-critical current values (J -C 1). The 
thermally assisted properties are then easily understood 
by studying the mean first passage time over the effec- 
tive barrier. This is done by contructing and solving the 
associated adjoint equation to the Fokker-Planck opera- 
tor [H, Q. In the limit of low noise (£ >■ 1) the scaling 
behavior between mean switching time and current for 
sub-critical currents becomes: 



y¥ F[ve(i-/)]cx P [e(i-/) 2 ] 

a I -I 2 



(4) 



where F[x] — exp(— x 2 ) exp(y 2 )rfy is Dawson's inte- 
gral. 

To verify this result, we numerically solved the full 
stochastic LLG equation (Eqn. 1) by using a standard 
second order Heun scheme to ensure proper convergence 
to the Stratonovich calculus As mentioned above, 

the nanomagnct is initially in thermal equilibrium at a 
temperature T in one of its two possible stable states at 
the start of the simulation (r = 0), when the current is 
applied. At each time step, the random kicks arc gener- 
ated with strength given by 



D 



ak B T 



2K(l + a 2 ) 2(l + a 2 )£ 



(5) 



which is obtained by imposing the fluctuation-dissipation 
theorem in the absence of applied currents. Such a diffu- 
sion constant is valid as long as the magnetic sample is 
in equilibrium with the thermal bath and effects due to 
Joule heating are neglected. 

Statistics were gathered from an ensemble of 5000 
events with a step size of 0.01 in reduced unit of time. 
The Landau damping constant is set at a = 0.04. Al- 
though the main results in this paper are shown for a 
barrier height of ^ = 80, different barrier heights were 
explored. To solve the numerical equations out to long 
time regimes, the event dynamics were simulated in par- 
allel on an NVidia Tesla C2050 graphics card. The large 
number of necessary random numbers were generated by 
employing a combination of the three-component com- 
bined Tausworthe "taus88" [TU and the 32-bit "Quick 
and Dirty" LCG [13] whose statistical properties are 
discussed in the literature This hybrid generator 

provides an overall period of 2 121 . Comparison of our 
result (4) with numerical data is shown in Fig. [T] a 
(t) cx exp(— £(1 — I) 2 ) is a much better fit to the data 
at long time scales than (r) cx cxp(— £(1 — /)). Nonethe- 
less, even at mean switching time scales four orders of 
magnitude larger than switching times at critical current, 
full asymptoticity in the numerics is still not reached 
exactly. This would seem to imply that the crossover 
regime, a regime characterized by an interplay between 
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FIG. 2: Comparison of the mean switching time versus cur- 
rent in the sub-critical low current regime for different val- 
ues of angular separation between polarizer and easy axes, u), 
ranging from to 7r/4. All critical currents have been rescaled 
to 1. Times are shown in units of (s • T) where T stands for 
Tesla: real time is obtained upon division by Hk- It is appar- 
ent that, up to a current rescaling, the mean switching time 
vs applied current is independent of uj. 



deterministic and thermal effects, spans over very large 
mean switching timcscales. 

We now consider how breaking the collinear symmetry 
impacts the scaling behavior. In particular, (3) no longer 
provides a simplification, because it depends explicitly on 
both m z and m y . These would need to be solved sepa- 
rately, given that q is now a third variable. Nonetheless, 
in the thermally dominated regime one can still make 
some useful approximations. For small values of a, the 
term in square brackets is of leading order over the second 
ballistic term depending on m y . This allows us, in the 
small a regime, to neglect the second ballistic term al- 
together. Furthermore, for low sub-critical currents, we 
expect an initially anti-parallel configuration to diffuse 
not very far from its local energy minima and establish a 
transient equilibrium until a sufficiently strong thermal 
kick drives it over the effective barrier over timcscales 
much smaller than those due to the deterministic spin- 
torque dynamics. This notion can be used to eliminate 
the m z dependence appearing in (3) by imposing its aver- 
age value over constant energy trajectories. These in turn 
are nearly circular orbits around the uniaxial anisotropy 
axis for which (m z ) = n z q. Doing this leads the second 
term in the square brackets to vanish until the onset of 
the switching event. We can then think of (3), in such a 
limit, to behave in the following approximate form: 



q = a(nj + q)(l - q 2 ) + J |(1 - q 2 ) o W. (6) 
This is reminiscent of the ID projcctional dynamics 
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discussed in relation to the collincar limit. The only dif- 
ference is the substitution I — > n z I. Furthermore, analo- 
gously to the collincar limit, one can again define a criti- 
cal current I c as the current above which all initial mag- 
netic configurations will switch deterministically. This 
has already been derived elsewhere in the literature Q 
and can be shown to give I c = l/cos(u>), where uj is the 
angular tilt between uniaxial and polarizer axes. This 
in turn automatically implies that n z I = I / I c . In other 
words, the thermally activated dynamics are the same 
for all angular separations up to a rescaling by their re- 
spective critical current. We then expect that the mean 
switching time dependences remain functionally identical 
to the collinear case for all uniaxial tilts. This has been 
confirmed by comparison with data from our simulations, 
and the results are shown in Fig. [5] It is of interest to 
point out how the effect of angular tilt in thermally ac- 
tivated spin-transfer reversal is quite different from what 
one gets in studying the similar phenomenon of field 
driven reversal. It has been shown, in fact, that upon 
introducing an angular separation between easy and ap- 
plied field axes, the scaling exponent chan ges abruptly in 
value from 2 for the collinear case to 3/2 [13 ]. 

These results have important implications for the anal- 
ysis of experimental data in which measurements of the 
switching time versus current pulse amplitude are used 
to determine the energy barrier to magnetization rever- 
sal. Clearly use of the correct asymptotic scaling form 
is essential to properly determine the energy barrier to 
reversal. The energy barrier, in turn, is very important 
in assessing the thermal stability of magnetic states of 
thin film elements that arc being developed for long term 
data storage in STT-MRAM. Further work should ad- 
dress how these results extend to systems with easy plane 
anisotropy and situations in which the nanomagnet has 
internal degrees of freedom, leading to a break down of 
the macrospin approximation. 

We also note that current flow is a source of shot noise, 



which at low frequencies acts like a white-noise source in 
much the same way as thermal noise. It is therefore in- 
teresting to understand when this additional source of 
noise plays a role. For a magnetic layer coupled to un- 
polarized leads, the current induced noise on the mag- 
netization dynamics was found to be m~p~7r^p V (20| . 
where V is the voltage drop across the magnetic layer, 
while Tl/Yr is a dimensionless ratio characterizing the 
coupling strength of the magnetic layer to the left (L) 
and right (R) leads. Thus the noise is maximal (V/4) 
for perfectly symmetrical couplings, and is smaller in the 
limit of highly asymmetric contacts. This basic behavior, 
and the order of magnitude of the effect, is not likely to 
be modified by polarized leads. We argue that the tem- 
peratures at which the experiments have been performed 
current noise effects are not important. The experiments 
arc performed at room temperature where T = 300 K. 
For an all metallic device, such as a spin- valve nanopillar, 
the couplings are nearly symmetrical and at the critical 
current a typical voltage drop across the magnetic layer 
is less than 10 /xV or, cquivalently, 1 K. For a magnetic 
tunnel junction device V can be ~ 1 V. However, in this 
case the coupling is asymmetric. One lead (L) forms 
a magnetic tunnel junction with the nanomagnet, while 
the other (R) a metallic contact. This gives Tr/Tl > 10 4 
and a relevant energy ~ 1 K, again far lower than room 
temperature. It appears that current induced noise can 
only be important at room temperature for a nanomag- 
net coupled symmetrically between two tunnel barriers. 
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